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A template bank to search for gravitational waves from inspiralling compact binaries 

I: physical models 
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Gravitational waves from coalescing compact binaries are searched using the matched filtering 
technique. As the model waveform depends on a number of parameters, it is necessary to filter the 
. . . data through a template bank covering the astrophysically interesting region of the parameter space. 

^^ ' The choice of templates is defined by the maximum allowed drop in signal-to-noise ratio due to the 

C^ ' discreteness of the template bank, fn this paper we describe the template-bank algorithm that was 

^^ , used in the analysis of data from the Laser Interferometer Gravitational Wave Observatory (LIGO) 

Cn ■ and GEO 600 detectors to search for signals from binaries consisting of non-spinning compact objects. 

Using Monte-Carlo simulations, we study the efBciency of the bank and show that its performance 
is satisfactory for the design sensitivity curves of ground-based interferometric gravitational wave 
detectors GEO 600, initial LIGO, advanced LIGO and Virgo. The bank is efficient to search for 
f^ ' various compact binaries such as binary primordial black holes, binary neutron stars, binary black 

Cn I holes, as well as a mixed binary consisting of a non-spinning black hole and a neutron star. 

lO ■ PACS numbers: 

> 

m 

o 

^ n 

^^ I Long baseline interferometrie gravitational wave detectors are either already taking sensitive data (LIGO [1[, 

^^O GEO 600 2] and TAMA jj|) or will be ready to do so in the near future (Virgo |j|). The most promising candi- 

^^ ' dates for these detectors are coalescing binaries consisting of compact objects such as neutron stars (NS) and/or black 
O I holes (BH) in close orbit, losing energy and angular momentum through the emission of gravitational radiation. 
O"^' In the past decade a lot of research has been carried out in modelling the late-time dynamics of compact binaries 

^ , and the emitted radiation. To this end post-Newtonian (hereafter PN) approximation has been used in which all 
(^ the relevant physical quantities are expressed as a perturbative series in the parameter v which is a measure of the 
relative velocity of the two component masses of the binary. The approximation is now available to a rather high 
order w^ in PN theory 5, 6, 7, 8, 9, lOj. These models are definitely accurate enough to serve as faithful templates for 
matched filtering the radiation from binary neutron stars (hereafter, BNS) and binaries consisting of primordial black 
holes (hereafter, PBHs) [ll|. This is because both BNS and binary PBHs merge outside the sensitive band of current 
ground-based instruments and strong relativistic effects that come into play at the time of merger but not properly 
modelled by the PN theory are unimportant for their observation. In the case of binary black holes (BBH), however, 
PN models are not quite accurate since these systems will be ultra-relativistic {v '^ 0.2-0.4) when their frequencies 
enter the band where the detectors are most sensitive and some of the BBH systems would even coalesce at frequencies 
where the detectors have the best sensitivity. There are now models that go beyond the PN theory to predict the 
waveform given out by these systems and the different models for the binary dynamics predict significantly different 
waveforms. 

A. Signal models used in the search 

Current searches for gravitational waves from inspiralling compact binaries employ template waveforms that are 
based on the adiabatic, non-adiabatic and phenomenological models. The adiabatic waveforms are obtained by using 
the PN expansions of the conserved energy E and flux ^ and solving the energy balance equation dE/dt = —T, 
which in turn leads to an evolution equation for the angular velocity a;(i) and hence the phasing of gravitational 
waves (/j(t) — (2'k)~^ J uj{t)dt. An implicit assumption made in writing the energy balance equation is that E does not 
change (appreciably) over an orbital timescale. There has been significant amount of activity to understand the poor 
behaviour [l3 . \13^ of the PN scheme resulting in new improved models that are based on the use of re-summation 
methods to accelerate the convergence of the PN expansions. The Pade re-summation method J^ uses rational 
polynomial approximations of the original Taylor expansions to improve the convergence of the PN series. The non- 
adiabatic models directly integrate the equations of motion (as opposed to using the energy balance equation) and 
there is no implicit conservation of energy that is used in the orbital dynamics as, for instance, in the effective-one- 
body (EOB) approach [Ij, [la, llfl Ii3 ■ The EOB maps the real two-body conservative dynamics onto an effective 
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one-body problem wherein a test mass moves in an effective background metric. 

The PN, Pade and EOB models produce physical waveforms, as opposed to the phenomenological waveforms 
[l8Lll9ll20| . The latter waveforms capture the main features of the physical waveforms but use unphysical parameters. 
Nevertheless the phenomenological waveforms themselves fit very well with both the adiabatic and non-adiabatic 
models, because the scheme uses greater number of free parameters and the parameters are extended to vary over a 
wider range than that which is physically allowed. 

B. Matched filtering and template banks 

Generally speaking, since the signal from a coalescing binary can be modelled reasonably well we can employ 
matched filtering, which is the optimal strategy to search for signals buried in stationary, Gaussian noise. However, 
real interferometer noise is neither stationary nor Gaussian, and we have to use additional consistency checks, such as 
the x^-test J2ll |22| , to discriminate true events from broadband transients of instrumental or environmental origin. 

The matched filtering performs a phase-coherent correlation of the data with a template which is a copy of the 
model waveform weighted down by the expected power spectral density of the data. The gravitational wave signal 
depends on a set of continuous parameters such as the amplitude, the masses mi, m2, and spins of the stars and a 
fiducial reference time (taken to be either "time-of- arrival" to or "epoch-of-coalescence" tc) and the corresponding 
orbital phase 0o- The parameters describing the inclination of the binary orbit and the location of the source on the 
sky are degenerate with other parameters, and are therefore, ignored in the search problem, although they can be 
deduced once the response of three or more non-co-located detectors to a gravitational wave event is known. The 
eccentricity may also be ignored, since it is expected to be negligible (due to orbital circularization) by the time the 
gravitational radiation enters the frequency band (10-1000 Hz) of a ground-based detector. In this paper we will only 
concentrate on non-spinning systems for which the dimension of the search is restricted to 4 parameters {mi, ?7i2,, tc 
and 00 ) excluding the additional 8 parameters to be used in a spinning system search (angular momentum of the 
system and the spin orientation of each object). 

Since the parameters of an incoming gravitational wave signal will not be known a priori, the data must be searched 
with a set of templates, called a template bank, designed to cover the whole parameter space. The distance between 
the templates in the parameter space is governed by the trade-off between computational power and loss in detection 
rate due to the discrete nature of the template bank. The spacing should be chosen so that the loss in signal-to-noise 
ratio due to the mis-match of the template with the signal does not deprecate the detectability of these sources. 

We employ a geometrical approach '23', '2^ |2^ |2a, [23| in constructing a template bank, namely we define a metric 
on the parameter space and use local flatness theorem to place templates at equal distances. The template bank could 
be seen naively as a uniform grid in the parameter space. However, to emphasize that this is not, in general a trivial 
task, we suggest to consider placing dots on a sphere at equal distances. Placing templates on the parameter space 
depends upon two crucial things: (i) geometrical properties of the signal manifold, whether it is flat or curved, and (ii) 
the coordinates chosen for template placement, which could be curvilinear even when the manifold is flat, or "almost 
flat" . Indeed, we need to choose a "Lorentzian-like" (or close to it) coordinate frame. The second point is similar to 
the problem of laying a regular grid on the plane using either Cartesian or polar coordinates. The situation is even 
more complicated because the metric, in general, depends on the waveform model and the PN order one chooses. 

C. Filtering with physical and phenomenological templates 

The different PN models differ signiflcantly from each other in their prediction of the phasing of the waves, making it 
necessary to search for binary black holes and neutron stars by using not one but all the different families of waveforms 
currently available. Luckily enough, as we will show in this paper and an accompanying paper |29 |. it turns out that 
we can use a single bank for all the different physical waveforms explored thus far. This bank is currently used by the 
GEO 600 and LIGO detectors to search for BNS, binary PBHs and BBH [HllSlp. The signal from a BNS sweeps 
through the detector bandwidth and the two neutron stars can even merge outside the sensitive bandwidth of the 
current ground-based detectors and the waveforms corresponding to different physical models agree reasonably well 
among each other (with overlaps very close to unity). For higher mass systems, e.g. BBH, the situation changes, the 
deviation between the various models getting more significant in the sensitive bandwidth of the detectors. 

Since for BBH the model waveforms from different families do not agree well with each other, a template bank that 
employs one PN family of templates will not be efficient to detect signals from another family. One solution is to filter 
the data through a number of banks with templates based on all the different models. If the reasonable assumption 
that at least one of the available PN models is close to the true gravitational wave signal bears out then our bank 
consisting of templates from all the different PN families should be efficient enough in capturing the coalescence 
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events. Unexpectedly, we have found that the bank works well as long as the template and signal are both from the 
same family meaning a single template bank might suffice to filter the data through different PN families. 

Another approach to this problem is to use a single template bank of phenomenological waveforms |lq . Buonanno, 
Chen and Vallisneri have designed phenomenological waveforms that are motivated by the structure of the PN signal 
in the Fourier domain which, for high mass systems, has good overlaps with all the known models. They are able 
to do this by (a) introducing (in the case of non-spinning binaries) two more parameters than are necessary in the 
case of physical waveforms and (b) allowing the parameters to vary over a wider, and in some cases unphysical, range 
than in the case of physical waveforms. The advantage in using the phenomenological templates is that they have 
been designed to give a good fit not only to known models but also to waveforms which could be "in between" these 
models. However, there is also a disadvantage in using these waveforms as they are very poor in estimating the 
physical parameters of the emitting system. Moreover, using two filters instead of one increases the computational 
cost and the threshold to be used. It produces more false alarm due to larger size of the corresponding template 
bank. So we might want to filter the data with candidate events through a bank of physical templates for the sake of 
parameter estimation. 

In a companion paper we shall define and study the efficiency of a phenomenological template bank to capture 
different physical models of gravitational wave signals using designed and real sensitivity curves of the various inter- 
ferometers. We have verified that such a template bank performs well with the minimal expected (or higher) efficiency 
for signals from BBH systems with total mass in the range 6-40 Mq for designed sensitivity. For sensitivity seen during 
the second Science run of the LIGO interferometers the bank covers the parameter space with good efficiency in the 
same range of the total mass and a lower frequency cutoff of 40 to 100 Hz. 

D. Organisation of the paper 

In the rest of this paper we will discuss only issues related to the physical template bank where both signals and 
templates are based on the physical models of the gravitational waveform from an inspiralling binary. In Sec. Hll 
and mil we briefly review, respectively, the Fourier-domain representation of the signal and the technique of matched 
filtering in a geometrical language, with the aim of introducing the notation and language used in Section IIVI to 
construct a template bank. Section [l VI constructs a generic template bank to search for binary inspirals from sources 
with a wide range of masses and applicable to any interferometer. The template bank specifically targets equal mass 
binaries but it is suitable for binaries with a small mass ratio as long as the precessional effects due to the spin of 
the larger body are unimportant. Section^discusses the results of Monte-Carlo simulations performed with different 
design sensitivity curves from a variety of interferometers. Section IVII summarizes the main results of this paper and 
discusses the application of the template bank introduced in this paper to capture other target models to be presented 
in a separate publication. 

II. INSPIRAL SIGNAL IN THE STATIONARY PHASE APPROXIMATION 

For the purpose of the template bank design it suffices to use a simple model of the signal. We shall use the Fourier 
representation of the standard PN waveform in the so-called restricted PN approximation. In this approximation, 
one neglects the PN corrections to the amplitude of the signal while the corrections to the phase are fully taken into 
account to the highest order possible/available. 

Let us begin with the time-domain representation of the waveform. The response of an interferometric gravitational 
wave detector to arbitrarily polarized waves from an inspiralling binary of total mass M = mi + m2, mass ratio 
1] = mim2 / M'^ , at a distance D is given by 

Ht) - ^^ [7TMf{t)f' cos[^(t) + ^c], (2.1) 

where it is assumed that the detector's motion is unimportant during the time when the signal sweeps across its 
bandwidth. Here, f(t) is the (invariant) instantaneous frequency of the signal measured by a remote observer, the 
phase of the signal (p{t) = 2tt J '^ f{t)dt is defined so that it is zero when the binary coalesces at time t = tc, fc is 
the phase of the signal at tc and A is a numerical constant whose value depends on the relative orientations of the 
interferometer and the binary orbit which when averaged over all angles is 2/5 17]. 

One can compute the Fourier transform h(f) of the waveform given in Eq. (|2.1|l using the stationary phase approx- 
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imation |ll25|: 



A M ^/^ I ^n r 

M/) = -7rT7T\7^r'/'e^P ^^if■,tc,Vc,Xk)+^^ (2.2) 



M'l/) = 27r/tc + ^c + 5I^fe/''""^'- (2-3) 

fc 

The parameters ipc and tc are the so-called extrinsic, or kinematic, parameters and they are defined by the relative 
orientation and position of the source and the detector and a fiducial time-of-coalescence of the binary. Intrinsic 
parameters A^ are the chirp parameters related to the two masses of the binary. At the 2PN order there are four 
non-zero chirp parameters, Ag, A2, A3, A4 (Ai = owing to the fact that there is no 0.5PN term in the PN expansion), 
given by 

3 , 5 /743 11 \ , -37ri/3 

Ao = ,oo„._..rN../3 ^ A2 = ———-— -H —r/ , A3 



1287/(7rAf )5/3 ' ~ 967rryM V336 4 7 ' ^~ StjIVP/^' 

15 /3058673 5429 617 2] 

^ ~ 64r;(7rAf)i/3 1^1016064 ^ lOOs'' '^ lu"^ J ' ^'' 

The chirp parameters are all functions of the two mass parameters - the total mass M and the symmetric mass ratio 
rj. Instead of the masses we can choose any two of the chirp parameters to characterize the signal. The choice (Aq, A3), 
in particular, is convenient since the masses can be computed explicitly in terms of these two parameters: 

Often, it is convenient to use chirptimes tq and T3 , in lieu of chirp parameters A^ , given by 



To 



(7^^f/i)-^/^ r3 = -^(7r^f/i)-2/^ (2.6) 



256TTfLv" ' 8/L77 



where /l is the lower-cutoff frequency at which the generation of the templates starts. The value assigned to /^ should 
be chosen carefully: while a low value increases the size of the template bank a high value decreases the signal-to-noise 
ratio. Section llV Bl discusses how t o op timally determine the lower frequency cutoff. Finally, the chirptimes can also 
be inverted in terms of the masses |33 : 

The correspondence between the two sets of parameters {mi , 7712) and (tq, t^) is illustrated in Fig. |21 

III. MATCHED FILTERING SEARCH FOR INSPIRALS 

A matched filter is an optimal linear filter to detect a signal of known shape in stationary, Gaussian noise. It can 
be introduced in several ways and here we shall follow the geometrical formalism 25, 26] of signal analysis which is 
most convenient for the construction of a template bank |2^. The geometrical picture is more easily motivated by 
considering the detector output as a discrete time-series rather than a continuous function. Without loss of generality 
we shall follow this standard procedure in summarizing the basic results from differential geometric approach to signal 
analysis but often we shall also write the equivalent formulas obtained in the continuum limit. Although the signal is 
considered as a discrete series in time, it will be treated as a continuous function of the signal parameters^, which we 
shall denote as ■&'^. The notation and language followed here closely resembles that in Refs. |26ll27l l3 ^ . 



^ The discussion in this Sec. is quite generic. Thus, we shall use the symbol 1? to denote all the parameters of the signal, both intrinsic 
and extrinsic. In specific applications it might be easier to treat the maximisation of the signal-to-noise ratio over some of the extrinsic 
parameters as we shall briefly discuss in the beginning of the next Section. 
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A. A Geometric approach to signal analysis 

The output x{t) of a detector is not recorded continuously: it is sampled discretely, usually at equal intervals Ai, 
chosen so as to obey the sampling theorem (see e.g., Ref. [SJI)- Thus, the detector output lasting for a time T and 
sampled at a rate fs (/s = (Ai)~^) is a time-series consisting of A^-samples, N = T/At. The set of all possible detector 
outputs forms an iV-dimensional vector space. A gravitational wave signal characterized by a set of p parameters •d'^ , 
/i — 0, l,...,p— 1, sampled in a similar way is also an A^-dimensional vector. However, the set of all signal vectors 
don't form a vector space, although they do form a manifold (a p-dimensional surface in an A^-dimensional vector 
space), with the parameters serving as a coordinate system. 

Further structure can be developed on this manifold by defining a metric. It turns out that the statistic of matched 
filtering naturally leads to the definition of a scalar product between vectors, which can then be used to induce a 
metric on the signal manifold. We shall now turn to a brief discussion of the matched filtering in geometrical language. 

Let Xk = rik + Sk [x{t) ~ nit) + s(t)] denote the output of a stationary^ detector consisting of instrumental and 
environmental noise background Uk [n{t)\ and a possible deterministic signal Sk [s{t)\. We treat Uk as a random 
time-series drawn from a large ensemble whose statistical properties are those of the detector noise. We use {uk) to 
denote the average over the ensemble and assume that n^ is a stationary Gaussian random process with zero mean. 
In the frequency domain we denote these quantities as : Xm — fim + Sm [x{f) — "-(/) + s(/)], where Fourier series Sm 
[s(/)] of the time series Sk [s(i)] is defined by 



N-l 
fe=0 



N-l „oo 

^ Sfce-2"""'^'/^, s(/) = / s{t)e-^'''f'dt. (3.1) 

7, n ■/— OO 



Next, let us introduce the Hermitian inner product of vectors. Given vectors g and h their inner product, denoted 
{g,h), is defined as 

^ ^hk Jo Oh[j) 

Here "*" denotes the complex conjugate and Sh is the one-sided noise power spectral density (PSD) given by 

{nkn*J = ^S^kSkm, {n{h)n*{f2)) = \Sh{h)5{h - h)- (3.3) 



The inner product naturally leads to the concept of the norm \\h\\ of a vector h : \\h\\ = \J{h, h). A vector h is said 
to be normalized if its norm is equal to unity. If the vector h is to begin with not normalized then we can define a 
vector with unit norm by 

(3.4) 



vwm 



It follows from the theory of hypothesis testing that in order to detect a known signal, but with an unknown set of 
parameters "i?^, in stationary Gaussian background noise n one must construct the statistic ^ given by 



^,,^^^^M^-)+^.^^W^(^^ ,(,,)) (3.5) 

and maximize it over the parameters "d^ . Here h is the template and q = h/ Sh is the optimal filter. However, the 
terms template and optimal filter are often used interchangeably to mean either of h or q. Writing the above equation 
in terms of q we see that the statistic ^ is the cross-correlation of the detector output with the optimal filter q and not 
just with the template h. Although at first this seems strange, upon closer inspection one can see why the optimal 
filter is not just a copy of the signal we are looking for but the signal divided by the noise PSD. The noise PSD weighs 



^ Stationarity is the statement that the statistical properties of a detector, i.e. the distribution of its output, remain the same as a function 
of time. Mathematically, one expresses stationarity by demanding that the correlation of the noise at two different instants depends 
only on the absolute time-difference of the two instants but not on the two instants. 
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the correlation differently at different frequencies, making greater contribution to the region of the spectrum where 
the sensitivity is larger (i.e., Sh is smaller) and smaller contribution where the sensitivity is weaker (i.e., Sh is larger). 
The signal-to-noise ratio (SNR) p is the mean of the statistic £_ divided by the square-root of its variance given by 

(0 ii^^h)) {.s,h) 



^/(F^^ ^ {{x,hf) - {{x,h)f ^(MO' 



where the last equality follows from Eq. (j^.^fl and our assumption that the noise is a stationary random process with 
zero mean. Thus, the effect of filtering a signal s with a template h is equivalent to projecting the signal s in the 
direction of /i, i.e. p — {s,h). This projection is, obviously, the greatest when h — s (the matched filter theorem), 
giving p = ■y/(s, s). When the template is not the same as the signal the SNR is less than the optimal value. That is, 
in general p < y^{s, s). 

B. The mis- match metric 

We can use the inner product Eq. (|3.2|) to induce a metric on the signal manifold. In the rest of this Section we 
shall assume that our templates are normalized, that is {h, h) = 1. The distance between two infinitesimally separated 
normalized templates on the signal manifold is given by 

||/i(i9'^-hd^^)-/i(^^)|P = IKd^f'lf 

= g^.d^f^dr, (3.7) 

where /i^ is the partial derivative of the signal h with respect to the parameter t?^ and Einstein's convention of 
summation over repeated indices is assumed. The quadratic form g^^ = (/i^, hn) defines the metric induced on the 
signal manifold. This metric is the same as the mismatch metric defined in Ref. |26Ll27| . Let 0{'d'^^'D^ -\- 5-3^) denote 
the overlap between two infinitesimally separated vectors h{i!}^) and h{d^ + (5i?^), that is 

0(1?'", z9^ + Si)") = {h{df -f Sd^"), /i(i?^)). (3.8) 

Since our templates are normalized we have 0{-d^,'d^) = 1. For vectors that are infinitesimally close to each other we 
can expand the overlap about its maximum value of 1 {S'd^ —>■ 0): 

0(1?^, ir + d-d") = 1 - A'VI^^ S-d^Sd" + ..., (3.9) 

where M^^. = -^ gjr^- The mismatch M between the two vectors is M =1- C)(t9^, ■(?'" + di?'^) = g u.u6-d^ S-d" ■ Thus, 
the metric on the manifold can be equivalently constructed using either of the following formulas ,35.] : 

5m'' = [^iiih^] (3.10) 

1 d^O 



2 d'di^d^'' ' 



(3.11) 

Next, we introduce the concept of the minimal match MM |23l l26L l27l| . In searching for a signal of known shape 
but unknown parameters one filters the data through a bank of templates. The templates in the bank are copies of 
the signal corresponding to a set of values z?f, i = 1,2 .... The parameters "d^ are chosen in such a manner that 
(a) the bank of templates has at least an overlap of MM with any signal whose parameters are in a given range 

minmax(/i(i?^),/i(?9f)) > MM, (3.12) 

and (b) the number of templates is the smallest for the chosen value of MM. The quantity MM is a very important 
notion in constructing a bank of templates. Indeed, if we assume that the incoming gravitational wave signal s(i?'^) 
is reproduced exactly by our model, so that s{'d^) — Ah{'d'^), then we could have a drop in the SNR due to the 
coarseness of the template bank. Indeed, if the template nearest to the signal is hiiS^) then the SNR is given by 

p = meix{s,h,) = 710(1?^, i9f) >Ax MM. (3.13) 

hi 
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In reality, however, the overlap between the template and the signal can be smaller than the minimal match since our 
template may not be an accurate representation of the signal. Recall that our templates are based on the solution to 
a PN expansion of the Einstein's equations. Thus, our template might not be sufficiently accurate in describing the 
fully general relativistic gravitational wave signal emitted by a binary during the final stages of the coalescence of two 
compact objects. A lot of effort in the last 15 years has been focused on improving the PN templates by (a) working 
out the expansion to higher orders, and (b) employing re-summation techniques that accelerate the convergence of 
the PN expansion. Current searches for gravitational waves deploy several families of template banks representing the 
different classes of the signal models. If the signal s and the template h are from different families then an important 
notion that will be useful in studying the efficiency of the template bank is the so-called fitting factor introduced in 
Ref. [33 and defined as: 

FF(s,/i(i?f)) = max(s,/i(z?f)) (3.14) 

The fitting factor incorporates the degree of faithfulness of our templates in detecting the expected signal s. 

Finally, given the minimal match and the metric on the parameter space one can estimate the number of templates 
required to cover the desired range of parameters using |26| | 



m[g,MM]^\2J ^ ™ j j V\9\dr, (3.15) 

where p is the number of independent parameters. 

Strictly speaking, thequadratic approximation used in Eg. 13.91 is valid for (5i?'' -^ |20| but for practical purposes 
the approximation is vafid for MM = 90% - 95%. 

C. Application to binary inspirals 

Let us now apply the geometrical language of the previous Section to the case of binary inspirals. Our aim is 
to derive an explicit expression for the metric and then use it to set up a bank of templates. The discussion in the 
previous Section illustrated the general method for obtaining the metric on a signal manifold whose dimension is equal 
to the number of parameters characterising the signal. It might at first, therefore, seem as though we would need 
templates in the full multi-dimensional signal parameter space. We need a bank of templates since we would not know 
before hand what the parameters of the signal are and have to construct the detection statistic for different possible 
values of the signal parameters, such that at least one template in our bank is close enough to any possible signal 
so as not to lose more than, say 5% of the SNR. However, templates are not always explicitly needed to maximize 
the SNR (as shown by Schutz |37|) in the case of the extrinsic parameters tc and ipc, thereby greatly reducing the 
problem of having to filter the data through a large template bank. Indeed, data analysts are always looking for ways 
to reduce the effective dimensionality of the template space. 

The maximization over the phase (pc is achieved via the decomposition of the signal into its quadratures: two 
templates with a 7r/2 phase difference, e.g. (pc = and (pc — 7r/2 J23ll37| 



p = max(a;, k{^c)) = J{x, h,{0))^ + {x, k{^/2)f. (3.16) 



max 

•PC PC 



Thus, we need only two templates to search in the (^-dimension. Similarly, the search for the fiducial time tc at which 
the binary merges can be efficiently found in practice via the fast Fourier transform 

maxp = max p(tc), p(ic) = ^ ^l^e^-"*^. (3.17) 

tc tc ^ bhm 

Thus, we need the metric and the templates only in the two-dimensional space of intrinsic parameters, which are taken 
to be the chirp parameters Tq and T3 instead of the masses mi and m2 of the binary. However, it is more convenient 
to begin with the metric in the three-dimensional space of {tc, tq, t^,) and then project out the coordinate tc- 
In practice, we can also work with the chirp parameters 61 and 62 : 

01 ^ 27r/iTo, 62 = 2^/iT3. (3.18) 
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We can use either of Eas. H3.10|) or H3.11|) to work out the metric. Owen |2g used Eq. (|3.11|) and maximized 
the metric over the phase ^pc to obtain an explicit expression for the metric ^ap in the three-dimensional space of 

{tc, 01,02) ■ 



7a/3 = 2 (^[V^"^/3] - ^[^^]J'i^p]) ' 



(3.19) 



where ^a is the derivative of the Fourier phase of the inspiral waveform with respect to the parameter 6", that is 
tpa = d'i'/dOa, and J is the moment functional [Sg of the noise PSD. For any function a{x) the moment functional 
J is defined as, 



1 r^" a{x)x~^/^ , 
where I{q) is the gth moment of the noise PSD defined by, 

j-xv -qjZ 
Jxi^ Sh(X) 



(3.20) 



(3.21) 



Here x = f / fo (also xu = /u//o Etnd xl = /l//o)i /o is a fiducial frequency chosen to control the range of numerical 
values of the functions entering the integrals, /l > is the lower-cutoff chosen so that the contribution to the integrals 
from frequencies smaller than /l is negligible and /u is the upper frequency cutoff corresponding to the system's last 
stable orbit frequency. The metric gmk on the subspacc of just the masses (equivalently, chirp times) is given by 
projecting out the coalescence time tc- 
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We shall now derive an explicit formula for the metric on the sub-space of masses by using the chirp parameters 9i 
and 02 from Eq. H3.18|l as our coordinates. 

The starting point of our derivation is the Fourier domain phase Eq. (|2.3|l . By writing the A^ 's in terms of Oi and 
O2 we get 
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where the constants a^m are given by: 

aoi 
617 o 



3 _ llTT _ 743 / 25 \ 

-, ^21-—, a22-^^l^l 



1/3 



031 



0-4,1 = -T-rT^ , 042 



5429 /257r 
5376 V~2~ 



1/3 



043 



3 

~2' 

^xl/3 

10838016 V 47r4 



15293365 



(3.24) 



The gradients -00 = d'^/dtc, V'm = d'^/d0,n enter the moment functionals which can be computed from the foregoing 
expression for the Fourier phase 



V'o = 27r/, ip„ 



fc=0 



ik'^ 



(fc-5)/3 



(3.25) 



where n is the PN order up to which the phase is known, or the PN order at which the metric is desired. Expansion 
coefficients ^mfe can be considered to be a (2 x n) matrix, which to 2PN order is given by: 



* = 



aoi 021/^2 + (022/3) (^2/^1)'^' a4i/(?i + a42/(3(0?02)'^')-(a43/3)(e2/ei)'/' 

. -021^1/^1 + (2a22/3) (^1/^2)'^' 031 -2041^1/^!- (042/3) (ei/0|)'^' + (4a43/3) (02/^1)'^' 



(3.26) 



Template placement ~ v. 1.4, 2006/07/18 16:01:49 Template placement for inspiral searches 

It is useful to note that the moment functional of polynomial functions a(x) = ^ a^x'^' are given by 



J 



^ttkx'' 



= ^afcJ(7-3fc), (3.27) 



where J{q) is the normalized moment given by J{q) = I{q)/I{7). Using the definition of the metric in Eq. H3.19|l and 
projecting out the parameter tc, we find 



ymi 



1 V^ * * i J(^7 J. -^ Tno i.wn9 -^ [J(9 - k) - J(4) J(12 - k)] [ J(9 - j) - J(4) J(12 - j)] \ 

(3.28) 

IV. TEMPLATE BANK BASED ON SPA MODEL 

In this Section we will discuss the problem of constructing a bank of templates based on the geometric formalism 
introduced earlier. To this end we shall use the specific signal model of an inspiral binary discussed in Sec. m This 
model uses a restricted PN approximation in which the PN amplitude corrections containing higher order harmonics 
are neglected. The end result of the construction of the template bank is a set of points in the parameter space of 
chirptimes or, equivalently, the masses. Each point is associated with a template built from a specific signal model. 

In principle, there is nothing in the construction of the bank that forbids us to use as our templates a model that is 
different from the one discussed in Sec. ^ Indeed, we shall show in a companion paper that although we have used a 
specific signal model in the construction of the bank, the same bank works for a wide variety of other signal models. 
In this paper we shall mainly focus on the SPA model, and only as an introduction to a more exhaustive study we 
shall present one case of a physical model based on the time domain Taylor model at 2PN, restricted to the BNS 
search. 

We validate the performance of our template bank and quantify its efficiency (see Sec. for a precise definition) 
using Monte Carlo simulations in which a number of signals, with their parameters chosen randomly, are generated 
and their best overlap with the bank of templates is computed. The algorithm described here is implemented in the 
LIGO Algorithms Library (LAL) 1391, a nd currently used by the LIGO Scientific Collaboration (LSC) to search for 
BNS, BBH and binary PBHs [lal^ls^. 

A. Effective dimensionality of the parameter space 

In order to construct the template bank we will start with the metric defined in Eq. (|3.7|l . As argued in Sec. IIIII 
the extrinsic parameters tc and Lpc do not need to be searched with a template bank. One can analytically maximize 
the overlap with respect to these parameters locally at each point in the space of masses. Thus, we project Lpc, tc 
onto the two-dimensional space of masses in which the chirptimes {tq,t^) defined by Eq. H2.6|l serve as coordinates. 
Instead of {tq^t^) we can equivalently take any two of the r^'s to be independent parameters to characterize the 
signal. The choice {tq,t^), is particularly attractive because their relationship to the masses is analytically invertible 
[331 [cf. Eqs . H2.5|l and 12.6|l ]. In addition, as we will argue later, the metric in these coordinates is a slowly varying 
function |41| . 

B. Lower frequency cutoff 

Before presenting an algorithm to place templates in the parameter space let us discuss how one can choose the 
lower frequency cutoff /l, which essentially determines the size of the parameter space of chirptimes and plays a 
crucial role in the computational resources required to process the data through the template bank. 

The initial fiducial frequency /^ defines the range of values of the chirptimes and is not itself a parameter to search 
for. However, it affects the length of the signals (therefore, the parameter space to be covered) and the SNR extracted. 

For example, the Newtonian chirptime tq oc /^ ; thus, lower values of /^ give longer templates with the immediate 
consequence of enhancing both the required storage space for templates and the computational cost to filter the data 
through the template bank. Furthermore, the number of cycles in a template increases as iVcyc oc [M-Jl)^^/^^ affecting 
both the overlap between different templates and the number of templates in the bank. 
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TABLE I: The number of templates required to search for BNS versus the lower cutoff frequency /l at a minimal match 
MM = 95%. As the cutoff is lowered the number of templates initially increases sharply but stabilizes after reaching a certain 
value depending on the noise characteristics of the detector in question. The last row gives the lower cutoff frequency that was 
chosen in our simulations by demanding that the loss in SNR due to this choice is not more than 1%. Note that the number 
of templates converges at smaller values of /l, but we don't gain in SNR by choosing /^ to be that at which the number of 
templates converge. 

/l(Hz) GEO 600 LIGQT Advanced LIGO Virgo 

10 6085 4382 18692 30650 

15 6099 4412 12782 22068 

20 6113 4449 9425 15992 

25 6129 4408 7086 12027 

30 6145 4312 5469 9140 

35 6002 3783 4312 6958 

40 4914 3121 3527 5504 

45 3768 2577 2885 4358 

50 3017 2271 2289 3458 

Chosen /l(Hz) 40 40 20 20 



What is really important, both for signal recovery and the number of templates required in the bank, is the effective 
number of cycles [131 which depends not only on the lower frequency cutoff but also on the noise PSD of the detector. 
The reason for this is the following: while the signal power increases with decreasing lower-cutoff as /^ , the noise 
PSD Shif) increases much faster below a certain frequency determined by the various noise backgrounds. Thus, there 
will be negligible contribution to the SNR integral and noise moments J{q) from frequencies below a certain lower 
cutoff. As a consequence, choosing the lower cutoff frequency to be smaller than a certain value has no advantage; 
on the contrary it has the undesirable effect of increasing the number of templates and the computational cost of the 
analysis. 

The constant minimal match MM has also to be set and we decided to chose a value MM — 95%. This value is 
large enough so that Eg. 13. 91 is still valid and small enough so that the number of templates fEa. l3.15ll is not too large 
with respect to computational cost. 

Table m summarizes the number of templates at a constant minimal match of MM — 95% as a function of /i. We 
can immediately see that for the different PSDs considered, the number of templates evolves significantly as Jl is 
reduced but stabilizes after reaching a certain value. Our aim is to choose the lower frequency cutoff small enough so 
that it has negligible effect on signal recovery. 

Choosing the lower frequency cutoff too high leads to a loss in the SNR because, as evident from Eq. (|2.2|l . there is 
greater power in the signal at lower frequencies, with the power spectrum of the signal falling with frequency as /^^' '^. 
However, the increase in signal power at lower frequencies does not help beyond a certain point because the SNR 
integrand is weighted down by the noise PSD of the instrument. Current generation of interferometers are dominated 
at frequencies below about 30-40 Hz by seismic noise. For instance, in the case of the Virgo interferometer the noise 
PSD raises as /"'*■* at frequencies below about 30Hz thereby dictating that setting /l ^ 30Hz might be a reasonable 
choice. 

Essentially, there are two opposing choices for /^ and one must make an arbitrary but optimal choice. In this paper, 
the lower frequency cutoff is chosen as large as possible with the constraint that no more than 1% of the overlap is 
lost as a result of choosing /^ to be different from 0. The choice of Jl would then depend, in principle, on the masses 
of the stars in the binary. This is because the templates are shutoff when the binary reaches the innermost circular 
orbit, making the bandwidth available to integrate the signal smaller, and forcing to choose a lower cutoff that is 
smaller for higher masses. However, again because of the steep increase in the seismic noise below a certain frequency, 
signals from binaries of total mass greater than ~ IOOM0 will not be visible in the first generation instruments. 

In Fig. n we plot the highest value of lower cutoff as a function of the total mass of the binary by requiring that the 
loss in the SNR is no more than 1%. The different curves in the plot correspond to the noise PSD expected in various 
ground-based interferometers. There is no need to make these plots for systems with different mass ratios since, at 
the level of approximation we are using, the cutoff depends only on the total mass and not on the mass ratio of the 
system. The left panel plots the lower cutoff frequency, which gives a loss in SNR of 1%, as a function of the total 
mass when the ending frequency of the system is set at the last stable orbit'^ (LSO). The right panel shows the same 



The last stable orbit frequeney is taken to be the one given for a test mass in Schwarzschild geometry: /lso = 1/(6"^' ^ttM) 
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result when the upper frequency is set at the light ring (EOB models at 2PN) . Our choice of lower frequency cutoff 
is summarized in Table Q] 

There is subtlety |43 in the choice of the lower-frequency cutoff. The optimal SNR is only sensitive to the lower- 
frequency cutoff but not the parameters of the template and the signal because they are both assumed to be the 
same. When filtering the data through a template bank, however, we have no optimal templates but for a set of 
measure zero signals. Consequently, what is important is not the behaviour of the SNR (which only depends on the 
moment J [7]) but that of the metric (which is sensitive to the difference in the parameter values of neighbouring 
templates). The metric is highly sensitive to the moment J[17], which converges rather slowly as a function of the 
lower-cutoff. We have chosen the easier approach described above in choosing the lower-cutoff, as that is less sensitive 
to the parameters of the system. The real justification comes from the fact that our choice of the lower-cutoff does 
not seem to have severely affected the efficiency of the template bank as we shall see later in this paper. 
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FIG. 1: The four curves in each plot show the lower frequency cutoff that corresponds to a loss of SNR of 1%, respectecly for 
the GEO 600, LIGO-I, Advanced LIGO and Virgo, using the ending frequency of the binary to be the last stable orbit (left 
panel) and the light ring (right panel). 



C. Template placement algorithm 



The templates are chosen in the parameter space such it is guaranteed that for every signal in the region of interest 
there will be at least one template in the bank with an overlap larger than or equal to the minimal match. In addition 
to the templates within the parameter space of the search a layer of templates might also be required just outside 
the boundary of the parameter space to assure the chosen minimal match for all signals within the parameter space. 
This happens because some of the templates that lie outside the boundary have coverage within the parameter space; 
without them signals with their parameters near the boundary might not have the desired minimal match '43| . 

The template placement algorithm consists of two steps: the first step deals with placing templates along a chosen 
curve in the parameter space, specifically along the equal-mass curve, while the second step deals with the placement 
of templates within and outside the boundary of the parameter space. 

The preference of chirptimes (or chirp parameters) over masses as coordinates on the signal manifold is dictated 
by the fact that these variables are almost Cartesian. This is clearly seen when considering the IPN model for which 
the metric g^i, = {h^, hi,) is a constant, i.e. independent of the chirptimes, while the same metric expressed in (Af, 77) 
or (mi,TO2) coordinates is not a constant. Thus, the signal manifold at IPN order is not only flat, the chirptimes 
are Cartesian-like coordinates. When we go to higher PN orders, the multi-dimensional signal manifold, in which all 
chirptimes are considered to be independent of one another, is again a flat manifold but the physical chirp manifold (i.e. 
the manifold formed by signals expected from the inspiral of black hole binaries) is only a two-dimensional manifold 
that is obtained by imposing the constraints T3 = T3(ro,Ti), T4 = T4(to,ti), and so on. This physical manifold, could, 
in principle, be a curved manifold but the curvature is not likely to be large since the absolute value and range of 
higher order chirptimes is small compared to the most dominant chirptimes. Hence the constraints are unlikely to 
result in a large curvature. We shall, therefore, assume that the metric is essentially constant in the local vicinity (on 
the scale defined by mismatch) of every point on the manifold. 
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FIG. 2: The left panel shows an example of the range of masses covered by the template bank while the right panel shows 
the corresponding region in the space of chirptimes tq-ts. The shape of the detector's PSD does not enter in the definition of 
the chirptimes but the lower frequency cutoff does and is taken to be /l — 20 Hz. The individual mass range is chosen to be 
[0.3, 30]Mo including the BNS, BBH and PBHs binary systems. 



In this work we shall specifically choose (tq, T3) to be the coordinate system. As we also mentioned above, the 
masses are analytically invertible in terms of this pair of chirptimes [cf. Eq. H2.5|) ]. The range of masses accessible 
to various interferometers is shown in Fig. |21 together with the corresponding space of chirptimes. Note that since 
rj < 1/4 not all of the M-77 or tq-ts space is physically meaningful. In the M-rj space the region above 77 = 1/4 is 
forbidden while in the space of T0-T3 the region below the curve marked mi — 7712 is forbidden"*. As we do not wish 
to place templates in this forbidden region (although this is, in principle, perfectly possible) the loss in SNR will 
be more than the minimal match if we do not place templates along the equal-mass curve. Imagine starting from 
the boundary on the left of the space of chirptimes and placing templates along the T3 — const, line. Eventually, 
the algorithm would take us beyond the boundary on the right and a proper coverage of the parameter space would 
require us to place a template in the forbidden region. Placing this last template on the equal-mass curve is not a 
solution as this could create a "hole" in the parameter space. For signals in the "hole" the overlap obtained will be 
less than the minimal match. 

Thus, our template placement algorithm consists of two steps: After computing the minimum and maximum chirp- 
times corresponding to the search space, namely (t™'", r™'''^) and (t™™, r™'*'^), the algorithm chooses a lattice of 
templates along the equal mass curve, lays a grid of templates in the rectangular region defined by the minimum and 
maximum values of the chirptimes and rejects the lattice point if the point itself lies outside the parameter space and 
none of the vertices of the rectangle inscribed within ambiguity ellipse lie within the parameter space. 



1. Templates along the 77 = l/A-curve 

In the first stage, templates are built along the equal mass (that is, rj = 1/4) curve starting from the minimum value 
of the Newtonian chirp-time and stopping at its maximum value. The algorithm is illustrated in Fig. |21 given the n-th 
template at O with parameters {tq,t^), and the distance {Stq,St^) between templates in our preferred coordinates, 
consider lines tq — Tq + Stq {QA in Fig. |3J) and T3 = t^ + 8t^ [PB in Fig. 13). The template next to {t^^t^), on the 
equal mass curve, must lie either along PB or along QA (cf. Fig.lJJ in order that all the signals that may lie on OAB 
are spanned by at least one of the two templates. Clearly, if we were to place the (n -I- l)-th template at B there will 
be a gap and some of the signals will not have the required minimal match; placing it at A meets our requirement. 

Note, however, that there is no guarantee that this will always work; it would indeed fail if the curve along which 



* Owen l43l has suggested that it might be worthwhile to express the signal entirely in terms of M and 77 (or entirely in terms of tq and 
T3) and explore the forbidden region. This is a way of expanding our net to catch the inspiral signals whose late time phase evolution 
is not acc urat ely described by the PN expansion. This suggestion was independently explored and extended by Buonanno, Chen and 
Vallisneri Il8i 
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templates are being laid is a rapidly varying function. However, there is no danger of this happening in the case of 
the 7] = 1/4 curve. Indeed, for increments Stq equal to the distance between templates expected at a minimal match 
of 0.95 we can assume the curve to be a straight line: T3 « Tj + t^Stq, where 7-3 is the derivative with respect to tq 
along the rj — 1/4 curve. Therefore, the above algorithm will never fail in our case. 
To locate the (n + l)-th template we compute the following pairs of chirptimes: 



n+l 






r- + Sr-, 
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11+ 1 



n+l 



4A, 



4Ao 



^n+l \ 2/5 



'0 



4Ao 

.n+l \ 5/2 



4A3 



(4.1) 



^0 



256(7r/i)8/3 ' ' 8(7r/L)5/3- 
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FIG. 3: Algorithm sketching the placement of templates along the 77 = 1/4 curve. The placement uses the rectangle inscribed 
in the ambiguity ellipse around each point of the parameter space considered. 



2. Templates in the rest of the parameter space 



In the second stage, the algorithm begins again at the point (t™'",t™'"), with the corresponding distance be- 
tween templates (Jtq"", (^r™'"), and chooses a rectangular lattice of templates in the region defined by (T™'",r™™) 
^^max^ .j-min-j j-^max^ ^max-j g^^^ j-^min^ ^max-j ^ j^ template is acccpted either if its coordinate is within the parameter 
space of search or if its span, for the given minimal match, has an overlap with the parameter space. It is this latter 
requirement that allows for a layer of templates just outside the boundary of the parameter space. These templates 
will have overlaps larger than or equal to MM for some signals within the parameter space of interest; leaving them 
out would cause 'holes' in the template bank where the match achieved by the template bank will be smaller than 
the minimal match. 

The implementation of the algorithm along the equal mass curve and in a rectangular lattice in the rest of the 
parameter space is plotted in Fig. ^ where the chosen templates are represented as points and the ellipses around the 
points. The ellipses define the regions where the overlap is greater than or equal to the minimal match. A sketch of 
the implementation of the algorithm is given in Appendix I VII for both the placement along the equal mass curve and 
the rest of the parameter space. 

In summary, the following parameters are required to compute the grid in the parameter space: 

• minimum mass of the component stars. 
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FIG. 4: Algorithm sketching the construction of a rectangular lattice of templates. Using the placement algorithm as shown in 
Fig. 1^ the parameter space can be spanned with ellipses using their inscribed rectangles. 



• maximum total mass of the binary system, 

• m,inimal m,atch which defines the distance between the nearby templates in the parameter space, 

• noise power spectral density S^if) of the detector in question, 

• lower and upper frequency cutoffs used in computing the moments of the noise PSD. 

V. VALIDATION OF THE TEMPLATE BANK 

We have carried out Monte-Carlo simulations to study the efficiency of the template bank described in the previous 
Section. By efficiency we mean the ability of the bank to capture signals with the loss in SNR no more than the one 
defined by MM, the minimal match. The Monte-Carlo simulations consist in generating a number Ns of random 
normalised signals s and finding the fitting factor FF between the signal and the bank of templates. We use s to 
denote the signal in order to distinguish it from the template h used in the bank. The models used to generate s and 
h are based on the SPA model at 2PN described in the Section |nj We randomize the parameters (7711,7712) in the 
range covered by the template bank as well as the initial phase (pc ■ Another issue is the randomization of the mass 
parameters, for which there are several choices: one choice could be to have a uniform distribution of component 
masses. However, this choice will not probe our template bank properly; indeed, the corresponding distribution of 
the total mass will be under-populated in the small or the large total mass range. In fact, in order to properly test a 
template bank, we should inject signals using uniform distribution of the parameters used to design the bank. Namely 
the set (to,T3). But again such a choice will over-populate the lower-end of the physical masses. We found it more 
convenient to inject signals with a uniform distribution in the total mass of the binary, although this choice is not 
necessarily astrophysical. Moreover, the total mass determines the frequency at which the PN waveform terminates 
and it would, therefore, be helpful in choosing our injections based on a uniform distribution in the total mass. Such 
a distribution would emphasize any problems with the efficiency of the template bank associated with waveforms that 
terminate in the sensitive bandwidth of the detector and/or have short duration. 

A. Figures of merit 

From the results of the Monte-Carlo simulations we construct several figures of merit (FOMs) that are relevant in 
validating the template bank. First, we define the efficiency of the template bank by a vector £ defined by 



£NAXs,Xh) ^ {FF {s,ir),hid'^))}^ 



^l..Ns 



(5.1) 
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where -d" and d^ are vectors corresponding to the parameters of the injected signals and the templates, Xs and 
Xh are the models used in the generation of the signal and template, respectively. In all the simulations, we set 
•d'^ = {toi, m,2, Veil but as mentioned earlier we can analytically maximise over the orbital phase tpc and, therefore, 
■ff'^ = {mi, 1712}. Moreover, we fix Xs = Xh = SPA model. 

From the efficiency vector £ and the injection parameter vector i)^, we can derive several FOMs 

1. £ versus the total mass. This is an important FOM which shows (1) if any injection has a value of FF less than 
the specified MM and (2) if so, for which range of total mass it occurs (see Fig. |HJ. 

2. A drawback of the first FOM is that it hides the information about the mass ratio 77 and its relationship to the 
efficiency. Although not important in the BNS case, where the minimum value of rj is equal to 0.1875, it can 
be interesting in other cases to look at the dependency of the efficiency on the 77 parameter. For instance, in 
the case of a binary consisting of a black hole and a neutron star (BH-NS) or BBH where rj can be significantly 
different from its maximum value of 1/4. 

3. A third interesting FOM is the cumulative distribution of £ as in Fig. [7| This is interesting since we can see 
how quickly the distribution drops when FF is close to the minimal match. If it drops sharply then it is a good 
sign telling us that the bank has the desired properties. However, a smooth drop is a warning of the coarseness 
of the bank or the presence of holes in it. Indeed the distribution is expected to be quadratic. 

Finally, let us introduce a quantity called safeness defined as 

S{Xs,Xh) =nim£{xs,Xh)- (5.2) 

Since we can perform only a finite number of simulations, we have to take into account the quantity Ng'. 

'5jVe (Xs,X/i) = min£:jv,(Xs,X/0- (5-3) 

This is the minimal FF found among the Ng injections. Higher the Ng the more confident we are with the value 
of the safeness. For instance, if Ng is set to 1 000, we might miss some holes in the template placement that might 
appear with a higher number of injections. If this quantity is less than the MM then the bank is under-efficient with 
respect to the minimal match chosen. Conversely, if the S ^ MM, then the bank is said to be efficient. The safeness 
is also dependant on the mass range considered. For instance, a bank can be efficient in a subset S of the mass range 
considered and under- efficient in its complement S . Any mass range in which the template bank is under-efficient 
should be considered separately and investigated more carefully. 



B. Simulation parameters 



The parameter of each simulations are 



• The mass range as defined by the minimum and maximum masses of the individual components of the binary. 
Although the parameter space as a whole need not be split into different categories, we found it convenient to 
study four different cases: BNS, BBH, BH-NS and binary PBHs, each having different number of templates and 
efficiencies. Our placement algorithm has been independently verified in the case of binary PBHs by D. Brown 
[Slllijl and found to be efficient. The mass ranges used in the different cases are 

— Binary PBHs: [O.S-ljAf©. In principle, we can choose the lower limit to be a value smaller than 0.3 M©, but 
in practice it implies a significant increase in the number of templates than we can handle. For example, in 
the case of LIGO-I, decreasing the lower limit to 0.2 Mq and 0.1 M© leads to roughly 130000 and 950 000 
templates, respectively. 

— BNS : [1-3] M©. 

— BBH : [3-30] M©. We choose the upper limit to be 30 Mq. The upper limit can be extended depending on 
/l; in the case of advanced LIGO or Virgo, we can expect to go up to 50 Mq. However, in our simulations 
we have chosen the same value for all the detectors. 

— BH-NS [1-30] Mq : we fix one of the component masses to be in the range [1-3] Mq and the other to be 
in the range [3-30] Mq. 
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TABLE IL The size of the template bank for the different searches and detectors. The BH-NS case includes here the BBH and 
BNS. See Fig. I^for a graphical representation. 

Bank size versus detector 

BBH 

BH-NS 

BNS 

Binary PBHs 

• Design sensitivity curves are GEO 600, LIGO-I, advanced LIGO and VIRGO as in Ref. [I3| have been used 
in our simulations. Since PSDs used are design sensitivity curves, the number of templates in each simulation 
remains the same as opposed to a Monte Carlo simulation involving real data (where the noise PSD would 
change from one data segment to the next hence changing the number of templates and their location in the 
parameter space). 

• The minimal match MM is set at 95% in all simulations. MM parameter directly affects the number of templates 
in the bank [cf . Eq. H3.15| . Table ^1 summarizes the number of templates in the bank for each detector and 
search. 

• For each combination of design sensitivity curve and type of search (BBH, BNS, . . . ) we carry out Ng = 10 000 
injections. 

• The sampling frequency is fixed at 4 kHz in all cases and, therefore, the templates and signals are shut-off either 
at the last stable orbit or the Nyquist frequency, whichever is smaller. 



C. Results and discussion 

Although the simulations we performed were split into 4 sub-categories, there is nothing which forbids us to cover 
the whole parameter space in a single category. However, in practice this has two major implications. First, for 
the full range of masses there are far too many templates. For instance, if we allow the template bank to cover the 
range defined by both binary PBHs and BBH then the number of templates become very large (160,000 in LIGO-I 
design sensitivity curve) . Second, splitting the template bank into sub-banks allows us to cover more extensively each 
specific binary and to tune the parameters used in each search. For instance, we can decrease the sampling frequency 
of the BBH search from 4 kHz to 2 kHz without any loss in the SNR because the signals in this region have an ending 
frequency below IkHz. 

It is not always intuitive to translate the parameter space into number of templates by looking only at the mass 
range. While scaling-laws for the number of templates in terms of chirp times and lower-frequency cutoff do exist 
[27||, the boundary effects are far too important when the parameter space volume is not too large compared to the 
volume coverage of each template. This is especially important towards the lower-end of the mass range where the 
two curves given by 77 = 1/4 and mi = m„iin (cf. dashed and solid curves in Fig.|21), meet and the template space is 
almost one-dimensional. For instance, extending the mass range of black holes from 30 Mq to say 40 Mq does not 
involve a significant increase in the number of templates. On the other hand, decreasing the lower mass of the binary 
PBH search has a deep implication: in the case of LIGO-I design sensitivity curve, the BBH mass range [3-30] Mq 
is covered with 774 templates; extending the range to 40 Mq changes the number of templates to 851. In the case 
of binary PBHs, going down from 0.3 Mq to 0.1 Mq increases the number of templates by a factor of 30. In Fig.[Sl 
we have shown the mass range for different astrophysical sources (LIGO PSD) and the number of template needed to 
have a 95% minimal match. The number of templates needed to cover the BBH area is of the order of a thousand. 
BNS and BHNS is four times more. And the PBHs case implies ten times more templates than the BNS case. 

We can extend the parameter to any possibility including hypothetical cases such as a binary composed of a PBH 
and a neutron star. The squares show binaries made of similar compact objects (BBH, BNS, binary PBHs) whereas 
the rectangles show cases where the component stars are different such as BH-NS, NS-PBH and BH-PBH binaries. 
The following simulations do not include the last two cases but we believe that our bank will work for these cases as 
well. 
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FIG. 5: The parameter space of the search can be split into sub-spaces corresponding to astrophysically interesting binary 
sources. The left panel shows the parameter space of masses of the component stars for various astrophysically interesting 
systems while the right panel shows the same in the space of chirptimes. The properties of the compact object and the number 
of templates needed to cover these sub-spaces are quite different. The solid triangles on the left plot show sub-spaces where 
both the objects are of the same type and in the same mass range. The dashed squares show binary systems in which the two 
objects belong to different ranges of mass and are possibly of different type. While binary PBHs, BNS, BBH and BH-NS are 
all astrophysically possible, binaries where one of the components is a PBH and the other is a neutron star or a black hole 
probably don't exist and need not be searched in the first instance. Note that in the space of chirptimes, where the area is 
proportional to the number of templates, these rather unrealistic systems occupy quite a large area and by excluding them in 
a search we can save significantly on the number of templates. 

1. BNS case 

Fh'st, we focus on the BNS case in which the mass range is set to 1-3 Mq. We plot bank efficiencies versus the 
total mass in Fig. El In these plots, we have combined the results of the 4 simulations involving GEO 600, LIGO-I, 
advanced LIGO and Virgo. In Fig. [7| we have separated the 4 results by plotting the cumulative histogram of the 
efficiencies. We can see that in all cases the safeness is above the minimal match: Si^i = 0.966, therefore the bank 
is efficient for the search of BNS. Furthermore, the cumulative histogram drops quickly while reaching the safeness. 
This is an indicator of the good behaviour of the bank. This also indicates that the bank is over-efficient and it is 
partly due to our choice of a rectangular, as opposed to a hexagonal, grid and partly due to the extra templates along 
the equal-mass curve. 

2. Binary PBHs 

The main difference compared to the BNS case is the number of templates. Indeed, we restrict the mass range to 
[0.3 — 1\Mq and, therefore, the number of templates increases by a factor of 10. Results are summarised in Fig. Eland 
Fig. [3 The conclusions are very similar to the BNS case. The safeness parameter is Si^i — 0.962 and therefore the 
bank is efficient for the search of binary PBHs. 



BBH case 



The BBH search requires far fewer templates than in the case of binary PBHs or BNS. However, it is important to 
recall that the duration of the signal decreases as the chirp mass increases and there is a critical mass above which 
no signal enters the detector bandwidth. In the following simulations we restrict the upper mass to be 30 Mq. 

By studying the results in Fig. Eland [3 and Table UTTl it is clear that the bank is efficient in the case of GEO 600, 
advanced LIGO and Virgo detectors, with the safeness equal to 0.973, 0.955 and 0.976, respectively, but in the LIGO-I 
case S < MM. From Fig. El we see that there are fewer than 1% of the injections that have a fitting factor below MM. 
However, these few injections do have an overlap larger than 0.930. The injections that fail to achieve the required 
match are all concentrated in a specific region of the parameter space with their total mass between 30 and 40 M© 
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FIG. 6: Distribution of the efficiency vectors versus total mass. The color bar on the right of each plot shows the number of 
recovered signals with given values of fitting factor and total mass. 

as seen from Fig. |S| This is the region where the approximation that the signal manifold is flat breaks down and 
causes a "small hole" in template placement; but the effect is insignificant and does not deserve a fix. In any case, it 
is possible to use a larger MM of the bank so as to capture these signals with the desired match. The increase in the 
number of templates is not a major issue since the number of templates in the BBH search is quite small compared 
to binary PBHs or BNS. 

4. BH-NS case 

Finally, we look at the BH-NS case. Here the component masses can take values in the range [1 — 30] M©. There 
is an important issue here related to both the injections and ranges of masses in the template bank. First, since the 
BNS and BBH cases have already been studied independently, we do not need to look again at BBH or BNS systems. 



TABLE III: Safeness parameter derived from the bank simulations. 
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FIG. 7: Panels from top-left to bottom-right summarize the cumulative distribution probability of the efficiency vector £ for 
the different ranges of masses considered (namely, BNS, BBH, BH-NS and binary PBHs). Each panel has results from all the 
different detectors. While BNS, BH-NS and binary PBHs show a quadratic drop when FF drops closer to the minimal match, 
the BBH case shows a slower drop in the case of LIGO-I and advanced LIGO, with some values below the minimal match. 

Therefore, we restrict nii G [1 — SJMq and TO2 G [3 — 60]Mq. Second, the bank has to be extended up to 63 Mq, 
which is the sum of the maximum mass of the two objects. The safeness in this case is above the minimal match with 
a value Si^i = 0.964, which means that the proposed bank is efficient for the BH-NS search as well. 

D. Taylor models 

Finally, we perform two simulations, in the BBH and BNS case using only LIGO-I design sensitivity curve, in which 
both the templates and signals are based on the standard time-domain post-Newtonian model at 2PN order jl3| . The 
parameters of the simulation are exactly the same as above. The results are shown in Fig. |3 As one can see, the 
recovered fitting factors are all above the MM for both BNS and BBH, which means that the bank designed in this 
paper using the SPA models can be used for physical models based on Taylor approximants as well. 



E. Computational cost 



The algorithm presented in this paper is relatively fast. The computation time required to generate a template 
bank is of the order of a second to a few tens of seconds on any standard computer, depending on the parameters 
considered. Table IIVI summarizes the computational cost needed to generate different template banks on a Linux 
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FIG. 8: The bank designed for the LIGO-I design sensitivity curve covers the mass range [3-30] Mq, but there is a 'hole' in the 
parameter space corresponding to the total mass in the range [30, 4:0]Mq. Black dots are where the templates are placed and 
crosses are the positions of the injected signals that have a fitting factor below the minimal match; majority of the injections 
are found with FF > MM and are not plotted for the sake of clarity. 

Atlhon 1.6 GHz machine. The parameters used in these examples are similar to those in Ref. [43|, that is Fl ~ 30 Hz, 
/., = 4096 Hz and MM = 0.95 and a VIRGO PSD. The Table shows that a BBH template bank is created in a 
second while a template bank with TTimin = O.SM© and JTimax — 30Mq (more than 180,000 templates) required only 
28 seconds. The number of templates quoted in Table Hvl and in Ref. 40] are in very good agreement if we take into 
account a 30% effect coming from the fact that we used a square lattice instead of an hexagonal one. 

Let us now consider our BBH template bank. The cost of matched filtering (to process 128 seconds of data through 
2422 templates) is of the order of 50 minutes while the template bank generation costs 1 second. It is obvious that the 
template bank generation is a small fraction of the total cost. Since the placement uses a square lattice the number 
of templates is expected to be an extra 30% more in comparison to an hexagonal grid and the computational cost is 
larger by the same factor. Nowadays, analysing a month of gravitational wave data through k, 2000 templates take a 
couple of days only depending on the number of computers used. 

Decreasing the number of templates is, therefore, interesting especially in the BNS and PBH cases where the number 
of templates is large although this is not a major issue in the BBH case where it might even be safer to use a square 
lattice (see Sec. IV(T3|) . We provide an improvement of the algorithm presented in this paper by incorporating an 
hexagonal placement in a future work |28|. 

VI. CONCLUSION 

In this paper we have used a geometrical method for constructing a template bank to search for gravitational waves 
from binary coalescences. The bank design is based on the Fourier-domain representation of the expected signal and 
the templates are placed on a square lattice. We have shown that the bank performs well with respect to its efficiency 
and safeness in the case of BNS, BH-NS and binary PBHs and for the 4 detectors considered and even over-efficient 
as expected since a square lattice has been used. 

The proposed bank is also efficient in the case of BBH when the total mass is in the range [6-30] Mq. It remains 
efficient above 30 Mq up to 60 in the case of VIRGO, advanced LIGO and GEO 600 detectors, the only problem 
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FIG. 9: Cumulative distribution probability of the efficiency vectors £ in the LIGO-I design sensitivity curve using a Taylor 
model instead of SPA both for templates and injected signals (top figures) and £ distribution versus the total mass (bottom). 

being the case of LIGO-I in a specific total mass range of [30-40] Mq. This particular case has been investigated 
and correspond to the largest masses considered where the length of the signal is very short. Increasing the minimal 
match from 0.95 to 0.97 in the area which gives low match cures the problem. Therefore, we conclude that the square 
bank proposed is efficient for all the PSDs and searches we studied in this paper. We remind here that this bank has 
been used by the LIGO Scientific Collaboration to search for binary PBHs and BNS in four science runs |29l ISOl Isif . 
Finally, it is important to emphasize that though we have been using signals and templates based on the SPA 
model, there is no restriction in using the same bank but for different models of signals and templates. Indeed, we 
have found that the current bank is efficient even when the signals and templates are both based on Taylor, Pade or 
EOB models (see Ref. [13 for model classification) as shown in a companion paper [23; in this paper we have also 
obtained a fast algorithm to construct a template bank on an hexagonal lattice with the number of templates reduced 
by « 30% (useful for the case of BNS and PBH searches where the number of templates could be very large). 
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APPENDIX A: PLACEMENT ALGORITHM 



In this appendix we provide the algorithm used for the template placement described in Sec. IIVI The algorithm is 
split into two distinct parts depending on the position in the parameter space. 

In a nut-shell, the algorithm to lay templates along the equal-mass curve is as follows: 

Begin at tq = t™''' 
do while (tq < t'^'''') 



T^ = ro + Sto, t^ = 4^3 (rcf /4Ao) 
ri = ra + 6T3, t^ = AAo (t3^/4A3) 
if ({t^,t^) is closer to (ro,T3) than {t^,t^)) 



f/5 
,5/2 
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{ 

} 
else 

{ 



To =T^,T3 =T^ 



To = To^,T3 ^ri^ 



} 

Add (tojTs) to InspiralTemplateList 

numTemplates++ 

Compute metric at (tojTs) 

Compute distance between templates at new point : ((5to, (Jts) 

} . 
The algorithm to lay templates in the rest of the parameter space is as follows: 



Begin at To = Tq^- , rg = Tg^'" 

Compute metric at (TojTa) 

Compute distance between templates at new point : ((5to, (Jts) 

Add (tojTs) to InspiralTemplateList 

numTemplates++ 

do while (T3 <= t^^"") 

{ 

Begin at tq = t™™ 

do while (to <= t^'^'') 

{ 

if ((tqjTs) is inside the parameter space) 

{ 

Compute metric at (tojTs) 

Compute distance between templates at new point : ((5ro, (Jts) 

Add (tqjT^) to InspiralTemplateList 

numTemplates++ 

} 

Increment tq : tq = tq + Stq 

} 

Increment T3 : T3 = T3 + 5t3 

Get next template along T3 = const. : (tq, r3) 

This algorithm is very simple to implement: except for the metric computation, the code is based on loops over 
simple calculations. It is a fast and robust algorithm amenable to easy implementation. 

ACKNOWLEDGMENTS 

This research was supported partly by Particle Physics and Astronomy Research Council, UK, grant PP/B500731. 
The authors benefitted from useful discussions with the members of the LIGO Scientific Collaboration, in particular 
B. Allen, P. Brady, D. Brown, J. Creighton and B.J. Owen. 



